Detecting the cosmic acceleration with current data 
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Abstract 

The deceleration parameter q as the diagnostic of the cosmological accelerating expansion 
is investigated. By expanding the luminosity distance to the fourth order of redshift and the 
so-called y-redshift in two redshift bins and fitting the SNIa data (Union2), the marginalized 
likelihood distribution of the current deceleration parameter shows that the cosmic acceleration 
is still increasing, but there might be a tendency that the cosmic acceleration will slow down in 
the near future. We also fit the Hubble evolution data together with SNIa data by expanding 
the Hubble parameter to the third order, showing that the present decelerating expansion is 
excluded within 2a error. Further exploration on this problem is also approached in a non- 
parametrization method by directly reconstructing the deceleration parameter from the distance 
modulus of SNIa, which depends neither on the validity of general relativity nor on the content of 
the universe or any assumption regarding cosmological parameters. More accurate observation 
datasets and more effective methods are still in need to make a clear answer on whether the 
cosmic acceleration will keep increasing or not. 
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I. INTRODUCTION 



More than ten years has passed since two independent groups discovered that our 
universe entered a stage of accelerating expansion using the Type la Supernova (SNIa) 
data pp. During the past years, many observations embrace the results, such as large scale 
structure [2], the cosmic microwave background (CMB) radiation [3J, and so on. All of 
these data strongly suggest that we live in an accelerating universe, with the deceleration 
parameter g < at low redshift, and an exotic component named dark energy is introduced 
to interpret this phenomena [1] . As more and more observational data released, the cosmic 
expansion history can be more accurately determined [SI E]. 

Besides the eternal acceleration scenario, such as the well-known ACDM model, some 
authors have proposed that a matter dominated decelerating expansion will resume soon 
after the acceleration starts, because the vacuum energy's anti-gravitational properties 
will reversed [7j. Using the CPL parametrization [8], Shafieloo et al. [9] recently found 
that the cosmic acceleration might be slowing down. This question is interesting, since it 
challenges the ACDM model, which fits most of the current data very well and predicts 
an eternal cosmic acceleration. But this result is far from final judgement. First, different 
parametrization methods of dark energy equation of state will show different influences 
on the evolution behavior of the universe [10] , and the CPL ansatz might be unable to fit 
the data simultaneously at low and high redshifts. Second, some tensions also might exist 
in the SNIa dataset which may lead to different results in the joint analysis [TT| 112]. with 
one suggesting a decreasing acceleration, others just the opposite. In order to avoid these 
disadvantages, one can consider the so-called cosmographic approach [T3|, by directly 
parameterizing the deceleration parameter, such as linear expansion q{z) = qo + qiz |14j . 
or other two-parameter model q{z) = | + [T5] . And some authors discovered that 

the universe might have already entered a decelerating expansion era by using only low 
(for example 2 < 0.1 or 2; < 0.2) redshift SNIa data [IHl IT7] . 

In this paper, we focus on the deceleration parameter as the diagnostic of the cosmo- 
logical accelerating expansion. Following the approach introduced in [16], we investigate 
the deceleration parameter by expanding the luminosity distance to the fourth order of 
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redshift and the so-called y-redshift [18], and fitting the latest SNIa sample (Union2) 
released by the Supernova Cosmology Project (SCP) Collaboration [19], which contains 
557 data including recent large samples from other surveys and uses SALT2 for SNIa 
light-curve fitting. But instead of cutting off the data at low redshift, we use the sub- 
samples with z < 1.0 throughout the paper, in order to circumvent the convergence issues 
discussed in [20] and to make our discussion physically meaningful, while applying some 
other technique to reduce the accumulation of systematic error, also we introduce some 
physical constraint conditions to make our constraint more robust [21]. To make an over- 
all comparison, the Hubble parameter H{z) is also expanded to the third order of the 
redshift and the ?/-redshift and fitted to the Hubble evolution data [221 123] • This consti- 
tutes the second part of the paper. In Sec. 3, we propose a non-parametrization method 
to reconstruct the deceleration parameters at differen redshifts directly from the SNIa 
data, and this method is independent of the calibration of the SNIa data. Sec. 4 contains 
our conclusion. 

II. COSMOGRAPHIC EVALUATION 
A. Methodology 

We adopt the assumption that the universe has a fiat Friedmann-Robertson- 
Walker(FRW) metric and then the scale factor can be approximated by the first several 
terms in a Taylor series; i.e., as a fifth order polynomial [23]. The coefficients in this 
expansion are considered as the parameters of the theory. Since one more term taken will 
greatly reduce the constraint ability of the parameter (see Table 1 in the first paper of [25] ) 
and four-parameter series can give a better approximation to the Hubble parameter and 
the distance modulus than other expansions [2T]. Therefore we take the first four terms 
only in our analysis. 

Accordingly, the luminosity distance ^^(z) can be expanded to the fourth order of 



3 



redshift [131 [El El], which is a generahzed form of the Hubble law, 







+ -(2 - 2go - 15go' - 15go' + lOgoJo + 5jo + So)z'] + 0{z'>) 



(1) 



where the coefficients are expressed with Ho,qo,jo, and Sq, the present values of Hubble, 
deceleration, jerk and snap parameters, respectively. This kinematic method can gives 
valuable information regarding the expansion rates of the universe, thus provides us with 
a testing ground for all cosmological solutions. To circumvent the convergence issues as 
discussed in [20], we choose the redshift cutoff z < 1.0 with 537 SNIa events involved. 

A so-called ^/-redshift is introduced in [20J, where y = In this case, the luminosity 
distance can be expanded as 



This expansion converges faster than the expansion with the redshift z and the system- 
atic error in this case is much smaller than the former case. This expansion has been 
extensively studied in [l6j . 




In principle, the luminosity distance (11 at low redshift {z < 1) can be exactly expressed 
by the infinite power series expansion of z or y, yet with infinite parameters involved. In 
practices, we choose the first four terms of the Taylor series as an approximation, which 
leads to the truncation error in our analysis. This is considered as the systematic error. 
The systematic error in ^ is relatively larger compared with the Union2 data. For 
example, when y = 1/2, the systematic error in [2] is about 3.1%, while it is 1.0% in the 
SNIa data. In order to make the systematic error smaller than that of the Union2 data, 
one has to cut off the data with z > 0.398. Here we employ another method to deal with 
this issue, rather than striking off high-;z SNIa data. First, we separate the data {z < 1.0) 
equally in two redshift bins and the divide is Zcut = 0.274, to guarantee that there are 
enough data in each bin. And then we employ Eqs. ([T| and ^ in each redshift bin, 
respectively. Note that for the z > Zcut case, the luminosity distance can be expanded 
around z = ^, which separates the data (0.274 < z < 1) equally in two parts. This 
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method will reduce the systematic error and give a better constraint. Accordingly, the 
luminosity distance for z > Zcut can be expressed as 

dL{z > Zcut) = ^[;^(2 - 2go - 15g^ - ISgjj + lOgoJo + Sjo + so){z - W 

-^(1 - go - + Jo - J(2 - 2go - 15g^ - ISgjj + lOgoJo + 5jo + so)) 

- + (^(2 - 2go - 15go' - 15go' + lOgoJo + 5jo + ^o) - ^ 
(1 - go - 3go' + Jo) + ^(1 - go))(^ - + (1 + ^(1 - go) 

-^(1 - go - 3go' + Jo) + ^(2 - 2go - ISg^ - 15go' + lOgojo + 5jo + so)) 

- I)] + ddl) ■ (3) 

Using the y-redshift, the divide is i/cut = 0.215 and the luminosity distance for the y > Ucut 
case can be expanded around the point y = \, 

ddy > Vent) = ^[^{50 - 26go + 21g2 - ISgjj + lOgojo - Tjo + so)iy - 
+ ^(^(50 - 26go + 21g2 - ISgj^ + lOgoJo - 7jo + So) 
+11 - 5go + 3ql - jo){y - + ^(^(11 - 5go + 3ql - jo) 
+ 1(50 - 26go + 21go' - ISgjj + lOgojo - 7jo + so) + 3 - qo){y - ^-f 
+(1 + ^(3 - go) + ^(11 - 5go + 3go' - jo) + ^ 
(50 - 26go + 21g2 - ISgj^ + lOgojo - 7jo + so){y - ^)] + d^i^) . (4) 

The systematic error of Eq. Q is about 0.411%, which is comparable with the most 
accurately detected SNIa data, whose error is 0.212%. 

The Hubble parameter H{z) can also be expressed with Ho,qo,jo,So by expanding it 
to the third order. Since q{z) is related to the second order derivative of dilz), and the 
first order derivative of H{z), this expansion is worthwhile to make a full analysis. Using 
the redshift and the y-redshift, the bubble parameter can be expanded as 

E{z) = ^ = 1 + (1 + go)^ + ^(jo - go')2;' + li^ql + H ' 3jo - 4gojo - so)^' , (5) 

tio I 
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and 

E{y) = ^ = l+(l+go)y+i(2+2go+jo-go)-2'+k6+6go-3go'+3go'+3jo-4goJo-so)2/'. 
Hq 2 

(6) 

Similarly, in order to reduce the systematic error, we expand the Hubble parameter in 
the regime of redshift z > z^ut as 

^^^^^ = [1 + go + ^(Jo - q'o) + ^(3go^ + 3ql - 3jo - 4goJo - ^o)](^ - ^) 
+ \[jo - go + ^(3go + 3go - 3jo - 4goJo - So)]{z - 
+ J(3go' + 3go' - 3jo - 4goJo - so){z - + ^ , (7) 

and 

= [1 + go + ^(2 + 2go + Jo - go) + ^(6 + 6go - 3ql + 3ql + 3jo - 4goJo 
-So)]{y - ^) + ^[2 + 2go + Jo - go + ^(6 + 6go - Sq^ + 3ql + 3jo - 4goJo 
-So)]iy - 1? + ^(6 + Qqo - 3ql + 3ql + 3jo - 4goJo - s,)){y - ^-f 



Ho 



(8) 



B. Datasets 



The dataset we use is the most recently released Union2 SNIa dataset [19], which 
contains 537 events for z < 1.0. We fit the SNIa by minimizing the value of the 
distance modulus and the xln Snia is obtained by comparing theoretical distance 
modulus fith{z) = 5\ogiQ[di{z)] + fiQ, where /iq = 42.384 — 51og;^Q h, with observed fiob of 
super novae: 

To reduce the effect of /io, we expand xL with respect to fiQ [28] : 

xL = A + 25/io + Ci^l (9) 
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where 



B 



A 




C 



1 




Eq. ([9]) has a minimum as 




sn,min 



A - B^C 



which is independent of /iq. In fact, it is equivalent to performing an uniform marginal- 
ization over /ig, the difference between xin ^^e marginahzed xL is just a constant 
|2H]- We will adopt xin goodness of fit between theoretical model and SNIa data. 

In order to discriminate the evolution behaviors of different models, another dataset is 
also employed, which is the recently summarized Hubble evolution data [22], containing 
7 data at redshift z < 1.0. We add in three more data a.t z = 0.24,0.34,0.43 by taking 
the BAO scale as a standard ruler in the radial direction [23], and also H{z = 0) = 
74.2 ± 3.6(A;m/(s ■ Mpc)) observed by the Hubble Space Telescope (HST) p6]. The Xh is 



To make use of the data, we perform a uniform marginalization over the nuisance 
parameter Hq and get 



defined as 
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a — c^/b, 



where 
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C. Results 



The analysis is performed by using the Monte Carlo Markov Chain in the multidimen- 
sional parameter space to derive the likelihood. It is natural to employ some physically 
obvious limitations to make the estimation of parameters more robust [2T]. We impose 
the positiveness conditions of the Hubble parameter and the luminosity distance when we 
move randomly in the parameter space, and select the points satisfying these conditions 
as our chains. Note that these conditions are independent of the content of the universe 
and the requirement of matter density fraction, In this sense they are robust and the 
priors are: 

d^z) > 0, 
H{z) > 0. 

We first investigate the constraint on the model parameters using the SNIa data, 
and add the Hubble evolution data for a comparison, then we analyze the evolution 
behavior of the deceleration parameters. The best-fitting values and errors of parameters 
for diiz), dL^y) and d)^ (z), d\ {y) (where the upper index H denote the case with 
Hubble evolution data added) are summarized in Table |l| in company with the best- 
fitting values of ACDM parameters. For the ACDM model, we use the same dataset, 
and the best-fitting value of current matter density fraction parameter Q^o = 0.268, thus, 
go = l^mo ~ 1 = —0.598. The marginalized likelihood distribution of go is also shown in 
Figure [T] and |2j respectively. 

If we use the SNIa data only, we find that the best-fitting values of the deceleration 
parameter are both negative, which are larger than that of the ACDM model, and we 
can not exclude the case of go > even within la error. What's more, the present value 
of the acceleration rate shows a slight tendency of decrease for the case of y-redshift 
expansion, compared with the case of redshift expansion. Note that even though the 
y-redshift expansion is systematically more accurate than the redshift expansion, but the 
figure of merit (FoM) of go, jo and so is smaller, it is 1.10 for redshift expansion and 0.02 
for y-redhift expansion. Compared to the redshift expansion, the y-redshift expansion 
gives a weaker constraint on the snap parameter. Recently, Xia et al. [23] investigated 
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the fifth order parameter cq using the similar cosmographic approach and found that the 
constraint on it is much weaker. 

With the Hubble evolution data added, the constraint ability is stronger, and the FoM 
for redshift and y-redshift expansion are 3.07 and 0.52, respectively. It worths noting that 
by employing the FoM to compare the constraint ability of the two methods, one often 
assumes that the parameters' log-likelihood surface can be quadratically approximated 
around the maximum, and the likelihood for the data is Gaussian. So, it may significantly 
mis-estimate the size of the errors if these assumptions do not hold. On the other hand, the 
Bayesian Evidence, E, provides us with a good method to compare different models |27j . 
which is not dependent on those assumptions. And the ratio of Evidences for two models 
Bi2 = E{Mi)/E{M2), also known as the Bayes factor, provides a measure with which 
to discriminate the models. We calculate the Bayesian Evidence ratio of the redshift 
and y-redshift expansion methods, and find that for the case of SNIa dada used only, 
InB = —0.83, while for the case with Hubble data added, InB = —0.37. It looks that the 

redshift expansion is a bit better. According to the Jeffreys grades, however, the result 
means that these two expansion methods are comparable in fitting the Union2 dataset. 
We see from the results that the redshift expansion approach excludes the case go > at 
about 95.4% confidence level, which tells us that the present universe is still in the stage 
of accelerating expansion. The constraint on the snap parameter is much tighter in the 
case of redshift expansion, and the constraint ability is improved in some degree with our 
method compared with some existing studies [131 ES] • 



It is of interest to investigate the evolution behavior of the deceleration parameter to 
compare with the ACDM model, since different behaviors may tell something on the 
universe. And another change in the sign of the cosmic acceleration is also interesting to 
us, if any. Using the best-fitting parameters, the deceleration parameter can be expanded 
with the redshift and y-redshift, respectively, as 

g(^) = go - (go + 2ql - jo)z + ^(2go + Sq^ + 8ql - 7goJo - 4jo - So)z\ (10) 
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TABLE I: The best-fitting values with \a and 2(t errors of c/O; Jo and sq. And the best-fitting 
values of ACDM model fitted by the same SNIa data. 




FIG. 1: ID marginalized distribution probability for the present deceleration parameter from 
the SNIa data. Curves are scaled to max[P{qo)] = 1. 

and 

Qiv) =qo- (go + 2go - jo)y + ^{Mo + ^Qo - ^^oJo - 2jo - So)y^. (11) 

We note that the quadratic form is essential for the detection of the transient acceleration 
phase, since in this scenario, q should be positive in the past {z > 2) and changes its 
sign at moderate redshift {zt ~ 0.3 — 1), then steps into the phase of g > 0. We plot the 
evolution behaviors of q compared with the ACDM model in Figure [3j Because of the 
weak constraint of the dL^y), the la area is too large to give any conclusive information, 
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FIG. 2: ID marginalized distribution probability for the present deceleration parameter from 
the SNIa data and the Hubble evolution data. Curves are scaled to max[P{qQ)] = 1. 

therefore we do not show it here. 



The best-fitting curves in the figures show a remarkably different evolution behaviors 
of the deceleration parameter from the ACDM model, all of which indicate a transient 
acceleration. The ACDM model predicts that the universe began to accelerate at z = 
0.35, but in our analysis, the redshift expansion predicts a earlier transition and the y- 
redshift expansion predicts the similar transition time as the ACDM model. However, 
all cases indicate a transient acceleration phase and predict a decelerated stage in the 
near future z < 0. This picture is even clear in the case of the y-redshift expansion 
if both datasets are used, which shows a considerably positive deceleration probability 
today and in the future, at 68.3% confidence level, and it indicates that the ACDM 
model is completely excluded in the future. The behavior of transient acceleration is also 
predicted or allowed by several dynamic models O |30] . Due to large uncertainties at high 
redshifts and in the future, however, it is still not able to confirm the existence of transient 
acceleration phase, nor the future decelerated acceleration, even at Icr confidence level. 

Next, we plot the curves of E{z) for the four cosmographic models and the standard 
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FIG. 3: Evolution behaviors of the deceleration parameter. The red one is the best-fitting curve 
of q{z) and the green line is the evolution of q{z) derived from the ACDM model. The blue 
area delimits the 68.3% confidence region for the q{z) reconstruction. 

ACDM model using the best-fitting parameters, compared with the Hubble evolution 
data in Figure |4j We see from the figure that the curves are all nicely consistent with 
the data, including the ACDM model. The two curves including the Hubble evolution 
data are much close to the one of the ACDM model than the other two. Note that 6 of 
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the 11 Hubble evolution data have large error bars, with cr(o6s) > 10%, and the worst 
data is a{obs) = 61.86%. So, better measured data are in urgent need to make a definite 
discrimination among different behaviors of our universe. 
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FIG. 4: Best-fitting curves of E(z) derived from the four cosmographic approaches. The red 
(or dashed red) line denotes the redshift (or y-redshift) expansion using SNIa data only. The 
green (or dashed green) line denotes the redshift (or y-redshift) expansion constrained with the 
Hubble data involved. The black one is the evolution of E{z) derived from the ACDM model. 



III. NON-PARAMETRIZATION APPROACH 

Further exploration on the deceleration parameter is also approached with a non- 
parametrization method, by reconstructing the deceleration parameter from the distance 
modulus of SNIa directly, which depends neither on the validity of general relativity nor 
the content of the universe or any assumption regarding cosmological parameters. From 
mathematics point of view, the determination of the deceleration parameter is much 
difficult than the evolution of Hubble parameter, since the Hubble parameter is the first 
order derivative of the comoving distance, while the deceleration parameter is related to 
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the second order derivative of the comoving distance. And because of that, the error for 
constructing the deceleration parameter is also larger. To deal with this issue, Ref. [6] 
generated 2000 artificial data sets mimicking what is expected from SNIa measurement by 
the SNAP satellite (see [22] or the web site ^) and then applied the quadratic expansion 
to express the dimensionless coordinate distance y{z) = HQr{z) in each redshift bin, thus 
the deceleration parameters can be obtained. 

There are two sources of errors we try to tackle in the implementation. First, the 
deviation from the true value of any individual data is assumed to be normally distributed, 
which is a good approximation only for relatively large samples. Second, any finite dataset, 
especially small dataset, will contain large sample variance. In summary, smaller number 
of data produces noisier results, while larger number of data leads to more convictive 
results at the expense of resolution in larger redshift bins. We try to balance these effects 
by choosing relatively large number of SNIa data in each redshift bin to make sure that 
there are both enough data and not too large bins, and also to reduce the effect of sample 
variance. 

In our analysis, we focus on the Union2 dataset instead of the mock datasets and 
separate the 557 SNIa data equally into four redshift bins, each containing 139 data 
points (149 data points in the forth bin). In each redshift bin, we reconstruct one data 
point using the data within each bin, so we get totally 4 reconstructed data, as shown 
in Table [IT], where the first column is the averaged redshift in each bin and the third 
column is the unbiased estimate of the error. We will not take use of the expression of the 
distance modulus of SNIa, which we believe contains some uncertainties in determining 
the calibration of the data. Instead, we apply the following expression, 

dL{Zi+i) 

/ii+i - /ii = 5 Ig — -, (12) 

where i runs from 1 to 3 in this case. This expression is credible if the distance modulus is 
of logarithm form of luminous distance only. In order to detect the deceleration parameter 
and make it easier for numerical approach, we expand the comoving distance to the second 

^ http://snap.lbl.gov 



14 



order of redshift in each redshift bin Zi < z < ^j+i, 



where the Hubble parameter can also be expressed with the deceleration parameter, 

= H{z,){\ + ^^(^.+1 - z,)\ (14) 
And then the luminosity distance can be obtained using the relation ^^(2;) = r(2;)(l + z), 

dL{z,^,) - - - ^(iT^^^^+i + TT^^^^^^^- ^^^^ 

For the beginning of the iteration process, we choose go ~ Q'l- In the following, we will see 
that this approximation is reasonable since the redshift distance between them is small, 
with 6z = 0.089. And we use the best-fitting value of Hubble parameter given by the 
WMAP7 group, ho = 0.704(/cm/(s ■ Mpc) [3\. Then we get 3 deceleration parameters. 



shown in Table |III[ where we choose the average redshift as the first column. Since 
q{z) involves the second order derivative of the distance modulus, it is much harder to 
make the constraint accurately, thus causes the error bars blowing up to be proportional to 
[211 [32], similar to the equation of state of dark energy w (2;). We plot the evolution 
of the reconstructed deceleration parameters in Figure [5j As a comparison, other four 
best-fitting curves of the cosmographic approach are also shown. It is clear that except 
the ACDM model, all other results contain a transient acceleration solution, though the 
present value of the deceleration parameter is still negative. We cannot discriminate their 
behaviors at low redshift {z < 0.4), but higher redshift behaviors are distinct {z > 0.6), 
thus more accurately detected high-z data are in urgent need. On the other hand, we still 
need a well behaved method to affirm if there is indeed a transient acceleration phase. 



IV. CONCLUSION 

To summarize, we have investigate the deceleration parameter in detail, both with the 
cosmographic approach and the non-parametrization method. These approaches depend 
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0.10286 


35.3573 


0.1867 


0.17760 


39.4708 


0.1616 


0.40451 


41.7502 


0.2926 


0.78596 


43.4330 


0.3439 



TABLE II: Reconstructed distance modulus of Union2 SNIa data, where the first column is the 
averaged redshift in each bin and the third column is the unbiased estimate of the error. 



z 




0-+ 


0"_ 


0.089 


-0.0834 


+2.0285 


-2.2051 


0.291 


-2.7717 


+3.1068 


-1.7184 


0.595 


2.2651 


+1.3969 


-0.9106 



TABLE III: Reconstructed deceleration parameters of Union2 SNIa data, where the last two 
columns are the estimates of the errors in each redshift bin, which are also directly constructed 
from the SNIa data by considering the full error of each data point. 

neither on the validity of general relativity nor the content of the universe. To make the 
results robust, we focus on reducing the systematic error throughout the paper. 

Cosmographic approach to detect the deceleration parameter is to expand the luminos- 
ity distance to the fourth order or even higher order of redshift or the so called y- redshift, 
and then to fit the datasets. But this method contains large uncertainties if we use the 
high- 2; data. Some authors cut off high-z data in data fitting process to avoid the problem, 
but we used other technique to reduce the systematic error: By expanding the luminosity 
distance in two redshift bins which contain the same number of SNIa data, the error can 
be reduced. Since the deceleration parameter is related to the second order derivative of 
the luminosity distance, and the first order derivative of the Hubble parameter, it is also 
necessary to expand the Hubble parameter to the third order as what we have done for 
the luminosity distance, thus improving the constraint ability. The result reveals that the 
universe may transit from decelerating expansion to a transient accelerating expansion 
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FIG. 5: Evolution of the deceleration parameter reconstructed from the Union2 SNIa data with 
la error bars. Others are the best-fitting curves of the cosmographic approach as before, in 
company with that of the ACDM model. 

during 0.3 < z < 1.0, and the best-fitting evolution behaviors tell us that it is possible 
that the universe will return to decelerating expansion stage in the future even though it 
is still accelerating at present. Comparing with the Hubble data shows that these results 
are indistinguishable with that of the ACDM model, which demands for more accurate 
observations to come up with a conclusive verdict. 

Non-parametrization method is also implemented by reconstructing the deceleration 
parameters directly from the SNIa data. We simply use the logarithmic form of the 
distance modulus with respect to the luminosity distance, thus avoid the uncertainties in 
calibration of the data. To reduce the systematic error, we distribute the same number 
of SNIa data in each redshift bin to make sure of enough data and not too large width 
of redshift bin. The result also indicates a transient accelerating expansion stage of our 
universe. Since the deceleration parameter is much more sensitive to the accuracy of the 
SNIa data than the Hubble parameter, it is not a easy job to constrain the result tightly, 
thus needs further improvement, both in observation data and fitting methods. 
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